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Abstract 

The relativistic viscous fluid equations describing the outflow of high tempera- 
ture matter created via Hawking radiation from microscopic black holes are solved 
numerically for a realistic equation of state. We focus on black holes with initial 
temperatures greater than 100 GeV and lifetimes less than 6 days. The spectra of 
direct photons and photons from decay are calculated for energies greater than 
1 GeV. We calculate the diffuse gamma ray spectrum from black holes distributed 
in our galactic halo. However, the most promising route for their observation is to 
search for point sources emitting gamma rays of ever-increasing energy. 

1 Introduction 



Hawking radiation from black holes fT]] is of fundamental interest because it relies on the 
application of relativistic quantum field theory in the presence of the strong field limit of 
gravity, a situation that could potentially be observed. It is also of great interest because of 
the temperatures involved. A black hole with mass M radiates thermally with a Hawking 
temperature = mp/SnM where mp = 

G-i/2 = 1,22 X 10^9 GeV is the Planck mass. 
(Units are h = c = k-^ = In order for the black hole to evaporate rather than accrete it 
must have a temperature greater than that of the present-day black-body radiation of the 
universe of 2.7 K = 2.3x 10~^ eV. This implies that M must be less than 1% of the mass of 
the Earth. Such small black holes most likely would have been formed primordially; there 
is no other mechanism known to form them. As the black hole radiates, its mass decreases 



1 



and its temperature increases until Th becomes comparable to the Planck mass, at which 
point the semi-classical calculation breaks down and the regime of full quantum gravity 
is entered. Only in two other situations are such enormous temperatures achievable: 
in the early universe and in central collisions of heavy nuclei like gold or lead. Even 
then only about T = 500 MeV is reached at the RHIC (Relativistic Heavy Ion Collider) 
just completed at Brookhaven National Laboratory and T = 1 GeV is expected at the 
LHC (Large Hadron Collider) at CERN to be completed in 2006. Supernovae and newly 
formed neutron stars only reach temperatures of a few tens of MeV. To set the scale from 
fundamental physics, we note that the spontaneously broken chiral symmetry of QCD gets 
restored in a phase transition/rapid crossover at a temperature around 170 MeV, while the 
spontaneously broken gauge symmetry in the electroweak sector of the standard model 
gets restored in a phase transition/rapid crossover at a temperature around 100 GeV. 
The fact that temperatures of the latter order of magnitude will never be achieved in 
a terrestrial experiment should motivate us to study the fate of microscopic black holes 
during the final days, hours and minutes of their lives when their temperatures have risen 
to 100 GeV and above. In this paper we shall focus on Hawking temperatures greater 
than 100 GeV. The fact that microscopic black holes have not yet been observed should 
not be viewed as a deterrent, but rather as a challenge for the new millennium! 

There is some uncertainty over whether the particles scatter from each other after 
being emitted, perhaps even enough to allow a fluid description of the wind coming from 
the black hole. Let us examine what might happen as the black hole mass decreases and 
the associated Hawking temperature increases. 

When Th <^ (electron mass) only photons, gravitons, and neutrinos will be created 
with any significant probability. These particles will not interact with each other but will 
be emitted into the surrounding space with the speed of light. Even when Th ~ me 
the Thomson cross section is too small to allow the photons to scatter very frequently 
in the rarified electron-positron plasma around the black hole. This may change when 
Th ~ 100 MeV when muons and charged pions are created in abundance. At somewhat 
higher temperatures hadrons are copiously produced and local thermal equilibrium may 
be achieved, although exactly how is an unsettled issue. Are hadrons emitted directly 
by the black hole? If so, they will be quite abundant at temperatures of order 150 MeV 
because their mass spectrum rises exponentially (Hagedorn growth as seen in the Particle 
Data Tables 0). Because they are so massive they move nonrelativistically and may form 
a very dense equihbrated gas around the black hole. But hadrons are composites of quarks 
and gluons, so perhaps quarks and gluon jets are emitted instead? These jets must decay 
into the observable hadrons on a typical proper length scale of 1 fm and a typical proper 
time scale of 1 fm/ c. This was first studied by MacGibbon and Webber 0] and MacGibbon 
and Carr p. Subsequently Heckler p argued that since the emitted quarks and gluons 
are so densely packed outside the event horizon they are not actually fragmenting into 
hadrons in vacuum but in something more like a quark-gluon plasma, so perhaps they 
thermalize. He also argued that QED bremsstrahlung and pair production were sufficient 
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to lead to a thermalized QED plasma when Th exceeded 45 GeV These results are 
somewhat controversial and need to be confirmed. The issue really is how to describe 
the emission of wavepackets via the Hawking mechanism when the emitted particles are 
(potentially) close enough to be mutually interacting. A more quantitative treatment of 
the particle interactions on a semiclassical level was carried out by Cline, Mostoslavsky 
and Servant They solved the relativistic Boltzmann equation with QCD and QED 
interactions in the relaxation-time approximation. It was found that significant particle 
scattering would lead to a photosphere though not perfect fluid flow. 

Rather than pursuing the Boltzmann transport equation one of us applied relativis- 
tic viscous fluid equations to the problem assuming sufficient particle interaction j^]. It 
was found that a self-consistent description emerges of a fluid just marginally kept in 
local thermal equilibrium, and that viscosity is a crucial element of the dynamics. The 
purpose of this paper is a more extensive analysis of these equations and their obser- 
vational consequences. The plan is as follows. In section 2 we give a brief review of 
Hawking radiation sufficient for our uses. In section 3 we give the set of relativistic 
viscous fluid equations necessary for this problem along with the assumptions that go 
into them. In section 4 we suggest a relatively simple parametrization of the equation 
of state for temperatures ranging from several MeV to well over 100 GeV. We also sug- 
gest a corresponding parametrization of the bulk and shear viscosites. In section 5 we 
solve the equations numerically, study the scaling behavior of the solutions, and check 
their physical self-consistency. In section 6 we estimate where the transition from viscous 
fluid flow to free-streaming takes place. In section 7 we calculate the instantaneous and 
time-integrated spectra of high energy photons from the two dominant sources: direct 
and neutral pion decay. In section 8 we study the diffuse gamma ray spectrum from 
microscopic black holes distributed in our galactic halo. We also study the systematics of 
gamma rays from an individual black hole, should we be so fortunate to observe one. We 
conclude the paper in section 9. 

2 Hawking Radiation 

There are at least two intuitive ways to think about Hawking radiation from black holes. 
One way is vacuum polarization. Particle-antiparticle pairs are continually popping in and 
out of the vacuum, usually with no observable effect. In the presence of matter, however, 
their effects can be observed. This is the origin of the Lamb effect first measured in 
atomic hydrogen in 1947. When pairs pop out of the vacuum near the event horizon of a 
black hole one of them may be captured by the black hole and the other by necessity of 
conservation laws will escape to infinity with positive energy. The black hole therefore has 
lost energy - it radiates. Due to the general principles of thermodynamics applied to black 
holes it is quite natural that it should radiate thermally. An intuitive argument that is 
more quantitative is based on the uncertainty principle. Suppose that we wish to confine 
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a massless particle to the vicinity of a black hole. Given that the average momentum of a 
massless particle at temperature T is approximately ttT, the uncertainty principle requires 
that confinement to a region the size of the Schwarzschild diameter places a restriction 
on the minimum value of the temperature. 

ttT ■ 2rs > 1/2 (1) 

The minimum is actually attained for the Hawking temperature. The various physical 
quantities are related as rs = 2M/mp = I/AhTh. 

The number of particles of spin s emitted with energy E per unit time is given by the 
formula 

dN, r, 1 

dEdt 271 exp{E/TH) - (-1)^' ' ^ ' 

All the computational effort really goes into calculating the absorption coefficient Tg from 
a relativistic wave equation in the presence of a black hole. Integrating over all particle 
species yields the luminosity. 



dt ' ' M2 

Here a{M) is a function reflecting the species of particles available for creation in the 
gravitational field of the black hole. It is generally sufficient to consider only those parti- 
cles with mass less than Th', more massive particles are exponentially suppressed by the 
Boltzmann factor. Then 



a 



2.011 X 10"^ [4200iVo + 2035iVi/2 + 835iVi + 95N2\ (4) 



where A^^ is the net number of polarization degrees of freedom for all particles with spin 
s and with mass less than Th- The coefficients for spin 1/2, 1 and 2 were computed 
by Page [|l^ and for spin by Sanchez |Tl|. In the standard model A^o = 4 (Higgs), 



A'i/2 = 90 (three generations of quarks and leptons), A^i = 24 (SU(3)xSU(2)xU(l) gauge 
theory), and N2 = 2 (gravitons). This assumes Th is greater than the temperature for the 
electroweak gauge symmetry restoration. Numerically a{TH > 100 GeV) = 4.43 x 10~^. 
Starting with a black hole of temperature Th, the time it takes to evaporate/explode is 

3a(r^)(87rT^)3- 

This is also the characteristic time scale for the rate of change of the luminosity of a black 
hole with temperature Th- 

At present a black hole will explode if Th > 2.7 K and correspondingly M < 4.6 x 10^^ 
g which is approximately 1% of the mass of the Earth. More massive black holes are cooler 
and therefore will absorb more matter and radiation than they radiate, hence grow with 
time. Taking into account emission of gravitons, photons, and neutrinos a critical mass 
black hole today has a Schwarszchild radius of 68 microns and a lifetime of 2 x 10^'^ years. 
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3 Relativistic Viscous Fluid Equations 



The relativistic imperfect fluid equations describing a steady-state, spherically symmetric 
flow with no net baryon number or electric charge and neglecting gravity (see below) are 
T^^^.j^ = black hole source. The nonvanishing components of the energy-momentum tensor 
in radial coordinates are [112 



J^00 
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7'(P + e)-P + t;'ATdi.s 
vYiP + e) + P + AT, 



diss 



(6) 



representing energy density, radial energy flux, and radial momentum flux, respectively, 
in the rest frame of the black hole. Here v is the radial velocity with 7 the corresponding 
Lorentz factor, u = vj, e and P are the local energy density and pressure, and 
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where rj is the shear viscosity and ( is the bulk viscosity. A thermodynamic identity gives 
Ts = P + e for zero chemical potentials, where T is temperature and s is entropy density. 
There are two independent differential equations of motion to solve for the functions T(r) 
and v{r). These may succinctly written as 



- (r'T^^) 
dr ^ ' 

— (r^T'-'-) 
dr ^ ' 




0. 



An integral form of these equations is sometimes more useful since it can readily 
incorporate the input luminosity from the black hole. The first represents the equality of 
the energy flux passing through a sphere of radius r with the luminosity of the black hole. 



(9) 



The second follows from integrating a linear combination of the differential equations. 
It represents the combined effects of the entropy from the black hole together with the 
increase of entropy due to viscosity. 
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The term L/Th arises from equating the entropy per unit time lost by the black hole 
—dS\^h/dt with that flowing into the matter. Using the area formula for the entropy of a 
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black hole, S^h = m|7rr| = AttM'^ /rrip, and identifying —dM/dt with the luminosity, the 
entropy input from the black hole is obtained. 

The above pair of equations are to be applied beginning at some radius tq greater than 
the Schwarzschild radius r^, that is, outside the quantum particle production region of the 
black hole. The radius tq at which the imperfect fluid equations are first applied should be 
chosen to be greater than the Schwarzschild radius, otherwise the computation of particle 
creation by the black hole would be invalid. It should not be too much greater, otherwise 
particle collisions would create more entropy than is accounted for by the equation above. 
The energy and entropy flux into the fluid come from quantum particle creation by the 
black hole at temperature Th- Gravitational effects are of order rs/r, hence negligible for 
r > (5 - lQ)rs. 



4 Equation of State and Transport Coefficients 



Determination of the equation of state as well as the two viscosities for temperatures 
ranging from MeV to TeV and more is a formidable task. Here we shall present some rel- 
atively simple parametrizations that seem to contain the essential physics. Improvements 
to these can certainly be made, but probably won't change the viscous fluid flow or the 
observational consequences very much. 

The hot shell of matter surrounding a primordial black hole provides a theoretical 
testing ground rivaled only by the big bang itself. To illustrate this we have plotted a 
semi-realistic parametrization of the equation of state in figure 1. Gravitons and neutrinos 
are not included. We assume (for fun) a second order electroweak phase transition at a 
temperature of Tew = 100 GeV. Above that temperature the standard model has 101.5 
effective massless bosonic degrees of freedom (as usual fermions count as 7/8 of a boson). 
We assume (also for fun) a first order QCD phase transition at a temperature of Tqcd 
= 170 MeV. The number of effective massless bosonic degrees of freedom changes from 
47.5 just above this critical temperature (u, d, s quarks and gluons) to 7.5 just below it 
(representing the effects of all the massive hadrons in the particle data tables) [|T^]. Below 
30 MeV only electrons, positrons, and photons remain, and finally below a few hundred 
keV only photons survive in any appreciable number. The explicit parametrization shown 
in figure 1 is as follows. 



A word about neutrinos: It is quite possible that they should be considered in approximate 
equilibrium at temperatures above 100 GeV where the electroweak symmetry is restored. 
Still there is some uncertainty about this. Since they provide only a few effective degrees 
of freedom out of more than 100 their neglect should cause negligible error. We will 
address neutrinos and their emission in a subsequent paper. 





Tew < T 

Tqcd <T< Tew 



QCD 



(11) 
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Now we turn to the viscosities. The shear viscosity was calculated in |H] for the full 



standard model in the symmetry restored phase, meaning temperatures above 100 GeV 
or so, using the relaxation time approximation. The result is 

77(T > 100 GeV) = 82.5 (12) 

when numerical values for coupling constants etc. are put in. The shear viscosity for 
QCD degrees of freedom only was calculated to leading order in the QCD coupling in 
^ to be 

^^^^^)-(l + iV,/6)a^ln(a-)^ ^''^ 

where Nf is the number of quark flavors whose mass is less than T. This is the extent 
of our knowledge of shear viscosity at high temperature. We observe that the ratio of 
the shear viscosity to the entropy density, as appropriate for the above two cases, is 
dimensionless and has about the same numerical value in both. Therefore, as a practical 
matter we assume that the shear viscosity always scales with the entropy density for all 
temperatures of interest. We take the constant of proportionality from the full standard 
model cited above. 

82.5 f s \ 

There is even less known about the bulk viscosity at the temperatures of interest to us. The 
bulk viscosity is zero for point particles with no internal degrees of freedom and with local 
interactions among them. In renormalizable quantum field theories the interactions are 
not strictly local. In particular, the coupling constants acquire temperature dependence 
according to the renormalization group. For example, to one loop order the QCD coupling 
has the functional dependence ~ 1/ ln(T/A) where A is the QCD scale. On account of 
this dependence the bulk viscosity is nonzero. We estimate that 

C^IQ-^V (15) 

and this is what we shall use in the numerics. 

Overall we have a modestly realistic description of the equation of state and the 
viscosities that are still a matter of theoretical uncertainty. One needs s{T),ri{T),({T) 
over a huge range of T. Of course, these are some of the quantities one hopes to obtain 
experimental information on from observations of exploding black holes. 



5 Numerical Solution and Scaling 

Several limiting cases of the relativistic viscous fluid equations were studied in The 
most realistic situation used the equation of state e = aT^, s = (4/3)aT^ and viscosities 
7] = hsT^i C = ^bT^ with the coefficients a, 65, hs all constant. A scaling solution, valid 
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at large radii when 7 1, was found to be T(r) = To(ro/r)^/^ and 7(r) = '^^{r /tqY^^ . 
The constants must be related by 36aToro = (32^5 + 4416^)70- This r-dependence of T 
and 7 is exactly what was conjectured by Heckler [0]. 

It was shown in |^ that if approximate local thermal equilibrium is achieved it can be 
maintained, at least for the semi-realistic situation described above. The requirement is 
that the inverse of the local volume expansion rate 9 = u^.^ be comparable to or greater 
than the relaxation time for thermal equilibrium . Expressed in terms of a local volume 
element V and proper time r it is ^ = {l/V)dV/dT, whereas in the rest frame of the black 
hole the same quantity can be expressed as {l/r'^)d{r'^u)/dr. Explicitly 

Of prime importance in achieving and maintaining local thermal equilibrium in a rela- 
tivistic plasma are multi-body processes such as 2 ^ 3 and 3^2, etc. This has been 
well-known when calculating quark-gluon plasma formation and evolution in high energy 
heavy ion collisions |16, 17| and has been emphasized in ref. [|], ^ in the context of black 



hole evaporation. This is a formidable task in the standard model with its 16 species of 
particles. Instead three estimates for the requirement that local thermal equilibrium be 
maintained were made. The first and simplest estimate is to require that the thermal 
DeBroglie wavelength of a massless particle, 1/3T, be less than 1/6. The second estimate 
is to require that the Debye screening length for each of the gauge groups in the standard 
model be less than 1/6. The Debye screening length is the inverse of the Debye screen- 
ing mass where n = 1,2,3 for the gauge groups U(l), SU(2), SU(3). Generically 
m° oc QnT where gn is the gauge coupling constant and the coefficient of proportionality 
is essentially the square root of the number of charge carriers |l^. For example, for color 



SU(3) = g^yl + N{/QT where A^f is the number of light quark flavors at the temper- 
ature T. The numerical values of the gauge couplings are: gi = 0.344, g2 = 0.637, and 
5^3 = 1.18 (evaluated at the scale mz) §]• So within a factor of about 2 we have mP ^ T. 
The third and most relevant estimate is the mean time between two-body collisions in 
the standard model for temperatures greater than the electroweak symmetry restoration 



temperature. This mean time was calculated in |TJ] in the process of calculating the 
viscosity in the relaxation time approximation. Averaged over all particle species in the 
standard model one may infer from that paper an average time of 3.7/T. Taking into 
account multi-body reactions would decrease that by about a factor of two to four. All 
three of these estimates are consistent within a factor of 2 or 3. The conclusion to be 
drawn is that local thermal equilibrium should be achieved when 6 < T. Once thermal 
equilibrium is achieved it is not lost because 6/T is independent of r. The picture that 
emerges is that of an imperfect fluid just marginally kept in local equilibrium by viscous 
forces. 

The results quoted above are only valid at large r and for the equation of state s ocT^. 
To know the behavior of the solution at non-asymptotic r and for the more sophisticated 
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equation of state and viscosities described in section 5 requires a numerical analysis. We 
have found that the most convenient form of the viscous fluid equations for numerical 
evaluation are 



47rr 

for energy conservation (from eq. (9)) and 



4 ( du u\ ( du 2u^ 

luTs - -ri-fu Ciu[ — + — 
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(17) 
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(18) 



for entropy flow (from eq. (10)). Obviously the entropy flux is a monotonically increasing 
function of r because of dissipation. 

Mathematically the above pair of equations apply for all r > 0, although physically 
we should only apply them beyond the Schwarzschild radius rs- Let us study them first 
in the limit r — > 0, which really means the assumption that v <^ 1. Then u v and 
7 f« 1. Wc also consider black hole temperatures greater than Tew so that the equation 
of state and the viscosities no longer change their functional forms. It is straightforward 
to check that a power solution satisfies the equations, with 

u{r) = Ui{r/rif'^ 

T{t) = T,{n/Tf"> (19) 

where is some reference radius. If the luminosity and the reference radius are given 
then Ui and Tj are determined by the fluid equations. 

The numerical solution for all radii needs some initial conditions. Typically we begin 
the solution at one-tenth the Schwarzschild radius. At this radius the Ui, as determined 
above, is small enough to serve as a good first estimate. However, it needs to be fine-tuned 
to give an acceptable solution at large r. For example, at large r there is an approximate 
but false solution: T — constant with u r-' r. The problem is that we need a solution 
valid over many orders of magnitude of r. If eq. (17) is divided by and if the right 
hand side is neglected in the limit r — oo then the left hand side is forced to be zero. 
We have used both Mathematica and a fourth order Runge-Kutta method with adaptive 
step-size to solve the equations. They give consistent results. 

In figure 2 we plot u{rs/ 10) versus Th- It is essentially constant for all Th > Tew 
with the value of 0.0415. In figure 3 we plot the function u{r) versus r for three different 
black hole temperatures. The radial variable r is expressed in units of its value when the 
temperature first reaches Tqcd, and u is expressed in units of its value at that same radius. 
This allows us to compare different black hole temperatures. To rather good accuracy 
these curves seem to be universal as they essentially lie on top of one another. The curves 
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are terminated when the temperature reaches 10 MeV. The function u{r) behaves hke 
r^^^ until temperatures of order 100 MeV are reached. The simple parametrization 

u(r)=us{r/rsf' (20) 

with Us = 0.10 will be very useful when studying radiation from the surface of the fluid. 

In figure 4 we plot the temperature in units of Tqcd versus the radius in units of 
tqcd for the same three black hole temperatures as in figure 3. Again the curves are 
terminated when the temperature drops to 10 MeV. The curves almost fall on top of 
one another but not perfectly. The temperature falls slightly slower than the power-law 
behavior r~^/^ expected on the basis of the equation of state s = (4/3)aT^. The reason 
is that the effective number of degrees of freedom is falling with the temperature. The 
entropy density is shown in figure 5. It also exhibits an imperfect degree of scaling similar 
to the temperature. 

Since viscosity plays such an important role in the outgoing fluid we should expect 
significant entropy production. In figure 6 we plot the entropy flow Anr US cLS cL function 
of radius for the same three black hole temperatures as in figures 3-5. It increases by 
several orders of magnitude. The fluid flow is far from isentropic. 



6 Onset of Free- Streaming 

Eventually the fluid expands so rapidly that the particles composing the fluid lose thermal 
contact with each other and begin free-streaming. In heavy ion physics this is referred 
to as thermal freeze-out, and in astronomy it is usually associated with the photosphere 
of a star. In the sections above we argued that thermal contact should occur for all 
particles, with the exception of gravitons and neutrinos, down to temperatures on the 
order of Tqcd- Below that temperature the arguments given no longer apply directly; for 
example, the relevant interactions are not those of perturbative QCD. 

Extensive studies have been made of the interactions among hadrons at finite tem- 



perature. Prakash et al. used experimental information to construct scattering am- 
plitudes for pions, kaons and nucleons and from them computed thermal relaxation rates. 
The relaxation time for vr — vr scattering can be read off from their figures and simply 
paramterized as 

/ T \^ 

^ 16 MeV (21) 

VlOOMeVy ^ ^ 

which is valid for 100 < T < 200 MeV. This rate is compared to the volume expansion 
rate 6 (see section 5) in figure 7. From the figure it is clear that pions cannot maintain 
thermal equilibrium much below 160 MeV or so. Since pions are the lightest hadrons and 
therefore the most abundant at low temperatures, it seems unlikely that other hadrons 
could maintain thermal equilibrium either. 
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Heckler has argued vigorously that electrons and photons should continue to interact 
down to temperatures on the order of the electron mass P, 0|. Multi-particle reactions 
are crucial to this analysis. Let us see how it applies to the present situation. Consider, 
for example, the cross section for ee ee'y. The energy-averaged cross section is ^ 

ai,,em = 8o;EMrl\n{2E/me) (22) 

where me is the electron mass, tq = aEM/^e is the classical electron radius, and E ^ nie 
is the energy of the incoming electrons in the center-of-momentum frame. (If one computes 
the rate for a photon produced with the specific energies O.IE, 0.25E, or 0.5E the cross 
section would be larger by a factor 4.73, 2.63, or 1.27, respectively.) The rate using the 
energy-averaged cross section is 



-1 

' brcm 



-e(3)T3 



ml 



ln(6T/me 



(23) 



where we have used the average energy {E) k, 3T for electrons with m^. -C T. This rate is 
also plotted in figure 7. It is large enough to maintain local thermal equilibrium down to 
temperatures on the order of 140 MeV. Of course, there are other electromagnetic many- 
particle reactions which would increase the overall rate. On the other hand, as pointed out 
by Heckler , these reactions are occurring in a high density plasma with the consequence 
that dispersion relations and interactions are renormalized by the medium. If one takes 
into account only renormalization of the electron mass, such that m^j^j ^ m^ + e^T^/3 
when me -C T, then the rate would be greatly reduced. 

Does this mean that photons and electrons are not in thermal equilibrium at the 
temperatures we have been discussing? Consider bremsstrahlung reactions in the QCD 
plasma. There are many 2 — > 3 reactions, such as: qiq2 — > (I1Q2I, <liQ2 ~^ 11921^ 99 ~^ QQl^ 
and so on. Here the subscripts label the quark flavor, which may or may not be the 
same. The rate for these can be estimated using known QED and QCD cross sections 



20| , pi| , P2| . Using an effective quark mass given by gT we find that the rate is a^T 
with a coefficient of order or larger than unity. Since becomes of order unity near 
Tqcd we conclude that photons are in equilibrium down to temperatures of that order at 
least. To make the matter even more complicated we must remember that the expansion 
rate 9 is based on a numerical solution of the viscous fluid equations which assume a 
constant proportionality between the shear and bulk viscosities and the entropy density. 
Although these proportionalities may be reasonable in QCD and electroweak plasmas at 
high temperatures they may fail at temperatures below Tqcd- The viscosities should be 
computed using the relaxation times for self-consistency of the transition from viscous 
fluid flow to free-streaming, which we have not done. For example, the first estimate for 
the shear viscosity for massless particles with short range interactions is T^r where r 
is the relaxation time. For pions we would get rj ~ constant, not rj ~ T^. As another 
example, we must realize that the bulk viscosity can become significant when the particles 
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can be excited internally. This is, in fact, the case for hadrons. Pions, kaons and nucleons 
are all the lowest mass hadrons each of which sits at the base of a tower of resonances HI • 



See, for example, and references therein. 

In order to do gamma ray phenomenology we need a practical criterion for the onset 
of free-streaming. We shall assume that this happens suddenly at a temperature Tf in the 
range 100 to 140 MeV. We shall assume that particles whose mass is significantly greater 
than Tf have all annihilated, leaving only photons, electrons, muons and pions. In figure 
8 we plot the freeze-out radius rj = r{Tf) for Tf = 100, 120 and 140 MeV versus the 
Schwarzschild radius. The fact that r/ increases as Tf decreases is an obvious consequence 

1/2 

of energy conservation. More interesting is the power-law scaling: rf^r^ ^ TjJ . This 
scaling can be understood as follows. 

The luminosity from the decoupling or freeze-out surface is 




(24) 



where the quantity in parentheses is the surface flux for one massless bosonic degree of 
freedom and df is the total number of effective massless bosonic degrees of freedom. For 
the particles listed above we have df = 12. By energy conservation this is to be equated 
with the Hawking formula for the black hole luminosity. 



L, = UT^^alffTl , (25) 

where a^^^ does not include the contribution from gravitons and neutrinos. Together with 
the scaling function for the flow velocity, eq. (20), we can solve for the radius 



and for the boost 



2 ( Ahnalff \"' Vt^ 
vr 1, 2u%df ) \ Tf 



(26) 



^^^^ = j " 0.22 . (27) 

From these we see that the final radius does indeed scale like the inverse of the square-root 
of the Schwarzschild radius or like the square-root of the black hole temperature, and that 
the average particle energy (proportional to 7/Ty) scales like the square-root of the black 
hole temperature. One important observational effect is that the average energy of the 
outgoing particles is reduced but their number is increased compared to direct Hawking 
emission into vacuum P, |^. 
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7 Photon Emission 



Photons emitted into the vacuum surrounding the black hole primarily come from one 
of two sources. Either they are emitted directly in the form of a boosted black-body 
spectrum, or they arise from neutral pion decay. We will consider each of these in turn. 



7.1 Direct photons 

Photons emitted directly have a Planck distribution in the local rest frame of the fluid. 
The phase space density is 

fiE') = ^E'/T, _ 1 ■ (28) 

The energy appearing here is related to the energy as measured in the rest frame of the 
black hole and to the angle of emission relative to the radial vector by 

E' ^-ff{l-VfCOS0)E . (29) 

No photons will emerge if the angle is greater than 7r/2. Therefore the instantaneous 
distribution is 



dEdt 



47rr5(— J d{cose) cos ef{E, cos 6) 



^ _M?>^ln(l-e-^/2^/^/) (30) 

where the second equahty holds in the limit 7/ S> 1. This limit is actually well satisfied 
for us and is used henceforth. 

The instantaneous spectrum can be integrated over the remaining lifetime of the black 
hole straightforwardly. The radius and boost are both known in terms of the Hawking 
temperature Th, and the time evolution of the latter is simply obtained from solving 
eq.(3). For a black hole that disappears at time t — we have 

Starting with a black hole whose temperature is Tq we obtain the spectrum 



dE TT^df \ 2uldf J E^ 



dxx^e""^. (32) 



Here we have ignored the small numerical difference between a^f^ and ah- In the high 
energy limit, namely, when E ^ 27j-(To)Tj, the summation yields the pure number 
4(27r^/315). Note the power- law behavior E'^. This has important observational conse- 
quences. 
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7.2 TT^ decay photons 

The neutral pion decays almost entirely into two photons: 7r° — > 77. In the rest frame of 
the pion the single photon Lorentz invariant distribution is 

^fN^ ^ 5{E-mJe) 

which is normalized to 2. This must be folded with the distribution of to obtain the 
total invariant photon distribution. 

d^pdt Jm.^ ^ dEj^dt TTrriT^ \ 2 / 

After integrating over angles we get 



rf^iV-" _ ^ i-oo dE^ d^N^, 



dEdt JE^in dE^^dt 



(35) 



where i?mm = {E'^ + ^"i/^)/ E. In the limit E ^ m^r we can approximate £^min — E and 
evaluate d'^N^o/dET^dt in the same way as photons. This leads to the relatively simple 
expression 

dEdt ~ ^.n'' ■ ^ ' 

The time-integrated spectrum is computed in the same way as for direct photons. 

dN^ _ 3604 /457r«^\ l_ r^/^vmr, 

dE - n^df [ 2uldf J E^ hn'h ^ ^ 

In the high energy limit, namely, when E ^ 2^f{To)Tf, the summation yields the pure 
number 27r^/315. 



7.3 Instantaneous and integrated photon spectra 

The instantaneous spectra of high energy gamma rays, arising from both direct emission 
and from 71^ decay, arc plotted in figures 9 (for Tf = 140 MeV) and 10 (for Tj = 100 
MeV). In each figure there are three curves corresponding to Hawking temperatures of 
100 GeV, 1 TeV and 10 TeV. The photon spectra are essentially exponential above a few 
GeV with inverse slope 2^f(TH)Tf oc ^TjTh- If these instantaneous spectra could be 
measured they would tell us a lot about the equation of state, the viscosities, and how 
energy is processed from first Hawking radiation to final observed gamma rays. Even the 
time evolution of the black hole luminosity and temperature could be inferred. 
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The time integrated spectra for Tj = 140 MeV are plotted in figure 11 for three initial 
temperatures Tq. A black hole with a Hawking temperature of 100 GeV has 5.4 days to 
live, a black hole with a Hawking temperature of 1 TeV has 7.7 minutes to live, and a 
black hole with a Hawking temperature of 10 TeV has only 1/2 second to live. The high 
energy gamma ray spectra are represented by 

dN _ mlTf 

dE ~ 26E^ • ^^^^ 

It is interesting that the contribution from tt^ decay comprises 20% of the total while direct 
photons contribute the remaining 80%. The fall-off is the same as that obtained by 
Heckler 0, whereas Halzen et al. p4| and MacGibbon and Carr |^ obtained an 



fall-off on the basis of direct fragmentation of quarks and gluons with no fluid flow and 
no photosphere. 



8 Observability of Gamma Rays 

The most obvious way to observe the explosion of a microscopic black hole is by high 
energy gamma rays. We consider their contribution to the diffuse gamma ray spectrum 
in subsection 8.1, and in subsection 8.2 we study the systematics of a single identifiable 
explosion. 



8.1 Diffuse spectra from the galactic halo 

Suppose that microscopic black holes were distributed about our galaxy in some fashion. 
Unless we were fortunate enough to be close to one so that we could observe its demise, 
we would have to rely on their contribution to the diffuse background spectrum of high 
energy gamma rays. 

The flux of photons with energy greater than 1 GeV at Earth can be computed from 
the results of section 7 together with the knowledge of the rate density p(x) of exploding 
black holes. It is 



dEdAdt 26 E^ J Aird^ix) 

where d{x.) is the distance from the black hole to the Earth. The exponential decay is due 
to absorption of the gamma ray by the black-body radiation The mean free path 
X^^{E) is highly energy dependent. It has a minimum of about 1 kpc around 1 PeV, and 
is greater than 10^ kpc for energies less than 100 TeV. 

We need a model for the rate density of exploding black holes. We shall assume they 
are distributed in the same way as the matter comprising the halo of our galaxy. Thus 
we take 
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where the galactic plane is the x — y plane, Rc is the core radius, and g is a flattening 
parameter. For numerical calculations we shall take the core radius to be 10 kpc. The 
Earth is located a distance Re = 8.5 kpc from the center of the galaxy and lies in the 
galactic plane. Therefore (f = {x — ReY + y"^ + z^- 

The last remaining quantity is the normalization of the rate density Pq. This is, of 
course, unknown since no one has ever knowingly observed a black hole explosion. The 
first observational limit was determined by Page and Hawking [^. They found that the 
local rate density piocai is less than 1 to 10 per cubic parsec per year on the basis of diffuse 
gamma rays with energies on the order of 100 MeV. This limit has not been lowered very 



much during the intervening twenty-five years. For example, Wright |23] used EGRET 
data to search for an anisotropic high-lattitude component of diffuse gamma rays in the 
energy range from 30 MeV to 100 GeV as a signal for steady emission of microscopic 
black holes. He concluded that piocai is less than about 0.4 per cubic parsec per year. 
(For an alternative point of view on the data see |28|.) In our numerical calculations 



we shall assume a value po = 1 pc yr corresponding to 

Plocal ~ 

0.58 pc ^ yr ^. 

This makes for easy scaling. Estimating the quantity of dark matter in our galaxy as 
-^haio/po, halo = 4.7 X 10^ kpc^ means that we could have up to 47 x 10^^ microscopic black 
hole explosions per year in our galaxy. 

Figure 12 shows the calculated flux at Earth, multiplied by E^. Of course this curve 
would be fiat if it weren't for absorption on the microwave background radiation. There 
is a relative suppression of three orders of magnitude centered between 10^^ and 10^^ eV. 
This means that it is unlikely to observe exploding black holes in the gamma ray spectrum 
above 10^^ eV. Even below that energy it is unlikely because they have not been observed 
at energies on the order of 100 MeV, and the spectrum falls faster than the primary cosmic 
ray spectrum oc i?^'^. The curve displayed in figure 12 assumes a spherical halo, q = 1, 
but there is hardly any difference when the halo is flattened to g = 2. 

8.2 Point source systematics 

Given the unfavorable situation for observing the effects of exploding microscopic black 
holes on the diffuse gamma ray spectrum, we now turn to the consequences for observing 
one directly. How far away could one be seen? Let us call that distance (imax- We assume 
that (imax < for simplicity, although that assumption can be relaxed if necessary. Let 
Adet denote the effective area of the detector that can measure gamma rays with energies 
equal to or greater than -Emin- The average number of gamma rays detected from a single 
explosion a distance dmax away is 



max " ^mm "^max 



Obviously we should have -Emin as small as possible to get the largest number, but it 
cannot be so small that the simple behavior of the emission spectrum is invalid. See 
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figure 11. 

A rough approximation to the number distribution of detected gamma rays is a Poisson 
distribution. 

P{N,) = ^^e-<^^) (42) 

The exact form of the number distribution is not so important. What is important is 
that when {N^{E > -Emin)) > 1 we should expect to see multiple gamma rays coming 
from the same point in the sky. Labehng these gamma rays according to the order in 
which they arrive, 1, 2, 3, etc. we would expect their energies to increase with time: 
El < E2 < < .... Such an observation would be remarkable, possibly unique, because 
astrophysical sources normally cool at late times. This would directly reflect the increasing 
Hawking temperature as the black hole explodes and disappears. 

It is interesting to know how the average gamma ray energy increases with time. Using 
eqs. (30) and (36) we compute the average energy of direct photons to be 47jTj(^(4)/^(3) 
and the average energy of tt" decay photons to be one-half that. The ratio of direct to 
decay photons turns out to be vr. Therefore the average gamma ray energy is 3.177/(t)Tj. 
This average is plotted in figure 13 for 10^ > t > 10~^ seconds. The average gamma ray 
energy ranges from about 4 to 160 GeV. 

The maximum distance can now be computed. Using some characteristic numbers we 

find 

i„ « 150^^ (-g— j PC. (43) 

If we take the local rate density of explosions to be 0.4 pc~'^ yr^^ then within 150 pc of 
Earth there would be 5 x 10^ explosions per year. These would be distributed isotropically 
in the sky. Still, it suggests that the direct observation of exploding black holes is feasible 
if they are near to the inferred upper limit to their abundance in our neighborhood. We 
should point out that a search for 1 s bursts of ultrahigh energy gamma rays from point 
sources by CYGNUS has placed an upper limit of 8.5 x 10^ pc~^ yr~^ [^. However, 
as we have seen in figure 13 and elsewhere, this is not what should be expected if our 
calculations bear any resemblance to reality. Rather than a burst, the luminosity and 
average gamma ray energy increase monotonically over a long period of time. 



9 Conclusion 

The increasing energy of the radiated photons by an exploding black hole and the dis- 
appearance of such a point gamma ray source in a certain period of time are unique 
characteristics of exploding black holes that may help us to detect them. Still, there is 
much work to be done in determining whether the matter surrounding a black hole can 
reach and maintain thermal equilibrium. The equation of state should be improved, and 
the viscosities computed using the relaxation times for self-consistency of the transition 
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from viscous fluid flow to free streaming. Also, there should be a more fundamental 
investigation of the relaxation times starting from the microscopic interactions. 

Our next step is to calculate the neutrino flux radiated by exploding black holes. 
Neutrino cross sections become very small at energies below 100 Gev which is the tem- 
perature of the electroweak phase transition. Above 100 GeV the neutrino cross sections 
are the same as other particles in the standard model which allows neutrinos to interact 
enough to reach thermal equilibrium. Therefore neutrinos are expected to freeze out at 
a temperature around 100 GeV. Another worthwhile project is to carry out cascade sim- 
ulations of the spherically expanding matter around the exploding black hole at a level 
of sophistication comparable to that of high energy heavy ion collisions. This project 
is much more complicated than the cascade simulation in heavy ion collision, though, 
because we need to deal with a much wider range of energies and particles involved in 
exploding black holes. 

The study of primordial black holes might well lead to great advancements in fun- 
damental physics. BecaTise the highest temperatures in the universe exist in primordial 
black holes, matter at extremely high temperatures can be studied, and physics beyond 
the standard model can be tested. In addition, because it is beheved that baryon num- 
ber is violated at high temperatures, the study of primordial black holes could possibly 
answer the question of why our universe became matter-dominated. Because primordial 
black holes explode, they are an ideal model for studying the Big Bang and the birth 
of our universe. Finally, the study of primordial black holes will help us to determine 
whether they are the source of the highest energy cosmic rays. The origin of these cosmic 
rays is still one of the biggest mysteries today. Observation and experimental detection 
of exploding black holes will be one of the great challenges in the new millennium. 
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Figure 1: Entropy density as a function of temperature, excluding neutrinos and gravitons. 
It is assumed that the QCD phase transition is first order and the EW phase transition 
is second order. 
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Figure 2: The value of u — at one-tenth the Schwarzschild radius as determined by 
numerical solution. The physical applicability of the numerical solution begins at radii 
greater than rs- 
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Figure 3: The radial dependence of u for three different Hawking temperatures. The 
curves begin at r^/lO and terminate when the local temperature reaches 10 MeV. 
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Figure 4: The radial dependence of T for three different Hawking temperatures. The 
curves begin at r^/lO and terminate when the local temperature reaches 10 MeV. 
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Figure 5: The radial dependence of s for three different Hawking temperatures. The 
curves begin at rs/lO and terminate when the local temperature reaches 10 MeV. 
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Figure 7: The rate for vrvr scattering and for the bremsstrahlung reaction ee — >■ ee7 are 
compared to the local volume expansion rate. The Hawking temperature is 10 TeV. 
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Figure 8: The freeze-out or free-streaming radius as a function of the Schwarzschild radius 
for three different freeze-out temperatures. 
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Figure 9: The instantaneous gamma ray spectrum for three different Hawking tempera- 
tures assuming Tf = 140 MeV. 
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Figure 10: Same as figure 9 but with Tf = 100 MeV. 
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Figure 11: The time-integrated gamma ray spectrum starting from the indicated Hawking 
temperature. Here Tf = 140 MeV. 
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Figure 12: The flux of diffuse gamma rays coming from our galactic halo. The normal- 
ization is /9o = 1 pc~^ yr^^. The halo is assumed to be spherically symmetric, q — 1] the 
results for a flattened halo with q = 2 are very similar. 
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Figure 13: The average gamma ray energy as a function of the remaining hfetime of the 
black hole. The times spanned correspond to approximately 400 GeV < Th < 200 TeV. 
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